Supernovae type la data favour coupled phantom energy 
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We estimate the constraints that the recent high-redshift sample of supernovae type la put on 
a phenomenological interaction between dark energy and dark matter. The interaction can be 
interpreted as arising from the time variation of the mass of dark matter particles. We find that 
the coupling correlates with the equation of state: roughly speaking, a negative coupling (in our 
sign convention) implies phantom energy (w<j, < — 1) while a positive coupling implies "ordinary" 
dark energy. The constraints from the current supernovae la Hubble diagram favour a negative 
coupling and an equation of state w$ < — 1. A zero or positive coupling is in fact unlikely at 99% 
c.l. (assuming constant equation of state); at the same time non-phantom values (w$ > — 1) are 
unlikely at 95%. We show also that the usual bounds on the energy density weaken considerably 
when the coupling is introduced: values as large as D. m o = 0.7 become acceptable for as concerns 
SNIa. We find that the rate of change of the mass m/m of the dark matter particles is constrained 
to be So in a Hubble time, with — 10 < <5q < — 1 to 95% C.L. We show that a large positive coupling 
might in principle avoid the future singularity known as "big rip" (occurring for w<j> < — 1) but the 
parameter region for this to occur is almost excluded by the data. We also forecast the constraints 
that can be obtained from future experiments, focusing on supernovae and baryon oscillations in 
the power spectra of deep redshift surveys. We show that the method of baryon oscillations holds 
the best potential to contrain the coupling. 



I. INTRODUCTION 

One of the main open problems in Cosmology is to 
determine the properties of the unclustered component 
called dark energy that is required to explain CMB, su- 
pernovae la, cluster masses and other observational data. 

Although most work on dark energy assumes it to be 
coupled to the other field only through gravity, there is by 
now a rich literature on possible interactions to standard 
fields and to dark matter. Even before the evidences in 
favor of dark energy, a coupling between a cosmic scalar 
field and ordinary matter has been studied by Wetterich 
P, Q while Damour, Gundlach and Gibbons Q inves- 
tigated the possibility of a species-dependent coupling 
within the context of Brans-Dicke models. Generally 
speaking, any low-energy limit of higher-dimensional the- 
ories predicts the existence of scalar fields coupled to mat- 
ter. More recently, several authors considered an explicit 
coupling between dark energy and matter: we denote this 
class of models as coupled dark energy. A partial list of 
works in this field is in refs. 0.0 |3| Other authors con- 
sidered a coupling to specific standard model fields: to 
the electromagnetic field Q , to neutrinos Q , to baryonic 
or leptonic current in these cases, there are strong 
constraints from local observations or from variation of 
fundamental constants. 

In this paper we confine our attention to the interac- 
tion with dark matter, as in ref. Q. Such a coupling 
is of course observable only with astrophysical 0] and 



cosmological ^lj experiments involving growth of pertur- 
bations or global geometric effects. Notwithstanding the 
large number of papers dealing with variants of coupled 
dark energy, the works dedicated to constraining the in- 
teraction via supernovae Type la (SNIa) are very limited. 
The main reason is that there is yet no compelling form of 
the interaction (just as there is no compelling form of the 
dark energy equation of state) . To make a step towards 
constraining the energy exchange between the dark com- 
ponents, instead of considering a specific coupling moti- 
vated by (or inspired by) some fundamental theory, we 
adopt here a phenomenological point of view. That is, 
we assume a general relation between dark energy and 
dark matter and derive its theoretical and observational 
properties. Our simple relation includes several previous 
coupled and uncoupled dark energy models but extends 
the analysis to cases which have not been tested so far. 
The main aim of this paper is to derive bounds on the 
interaction strength by analysing the recent SNIa sample 
of Riess et al. [13, which include several z > 1 super- 
novae. This works extends previous investigations in ref. 
[T^ |. which studied the same relation between dark mat- 
ter and dark energy, and in ref. > which compared the 
recent SNIa data with a special class of coupled dark en- 
ergy models motivated by superstrings and characterized 
by a constant ratio of densities. 



II. MODELING THE INTERACTION 
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Let us start with a model containing only dark matter 
and dark energy. The basic assumption of this paper is 
that the dark components interact through an energy ex- 
change term. The conservation equations in FRW metric 



2 



with scale factor a can be written in all generality as 



Pm + 3Hp m = SHp m , 



(1) 



where the subscript m stands for dark matter while the 
subscript <j> stands for dark energy and where 5 is a di- 
mensionless coupling function. In principle the coupling 
5 might depend on all degrees of freedom of the two com- 
ponents. However, if S is a function of the scale factor a 
only then the first equation of Q can be integrated out 
to give (we put the present value Qq = 1) 



Pm = Pmoa e 



J Sdc 



(2) 



where a = log a. This shows that the interaction causes 
p m to deviate from the standard scaling a -3 . That is, 
matter is not "conserved" or, equivalently, the mass m of 
matter particles p m = m(a)n, where n is their number 
density, varies with time such that 



m 
m 



(3) 



where the prime denotes derivation with respect to a. 
Therefore, S can be interpreted as the rate of change of 
the particle mass per Hubble time. 

While so far many paper have been devoted to find 
constraints on w^(a), aim of this paper is to derive con- 
straints on <5(a) from SNIa. We need therefore a sen- 
sible parametrization of 6(a). Generally speaking, we 
are faced with two possibilities: either we write down a 
parametrization for 6(a) and then find from this the rela- 
tion between p m and p^, or first give the latter and derive 
a function 5(a). We find that this second choice is simpler 
and better connected to previous work and constraints. 

The basic relation we start from is that the fluid den- 
sities scale according to the following relation 



P4>/p m = Acfi 



(4) 



where A,£ are two constant parameters. Moreover, we 
approximate as a constant: it is clear however that a 
more complete analysis should allow for a time-dependent 
equation of state. The relation J3J has two useful proper- 
ties: a) it includes all scaling solutions (defined as those 
with p ~ a m ) and 6) the functions p m (a), p$(a) can be 
calculated analytically. Since for uncoupled dark energy 
models with constant equation of state one has p m ~ a~ 3 
and P/j, ~ a^^ 1+w ^ it appears that the relation J3J re- 
duces to this case for £ = — iw<f>. Conversely, if £ — 3w^ 
then the matter density deviates from the a -3 law, as we 
have seen above. 

The relation Q has been first introduced and tested 
in ref. 0|. The present paper extends their investiga- 
tion in several respects. First, we include baryons. Ref. 
[l3T | assumed a single matter component, but there are 
strong upper limits to an interaction of dark energy with 
baryons (i.e. upper limits to a non-conservation relation 
like Eq. (2J for baryons, see e.g. flfilp. We will assume 



therefore that baryons are uncoupled while the coupling 
to dark matter is left as a free parameter. This has impor- 
tant consequences for the general behavior of the model. 
Second, ref. confined the range of parameters to 

< £ < 3 and w^, > —1, while we extend the range to a 
much larger domain, thereby including also the regime of 
"phantom" dark energy [JJ. Third, we use the new data 
of Riess et al. ^3|: these include SN at z > 1 and ex- 
tend therefore by a considerable factor the leverage arm 
of the method. Fourth, we'll discuss the implication of 
the results for as concerns the beginning of the acceler- 
ation. In ^3] it was claimed that the new high-redshift 
SN constrain z acc < 1 (more exactly, z acc = 0.46 ±0.13). 
In contrast with this, we show that z acc > 1 is not ruled 
out at more than 95% when the coupling is non-zero. 
A similar conclusion is obtained in ref. [lj| studying a 
string-motivated form of interaction and in adopt- 
ing different parametrizations of the equation of state. 
Finally, we will produce forecasts for future experiments. 

Before we include the baryons, let us derive some re- 
lations in the simplified case of two coupled components 
only. Assuming Sl tot = SI* + Sl m = 1 we obtain from J3J : 
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so that the constant A can be expressed in function of 
Sl^fi (the subscript indicates the present epoch) 
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n 
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If the dark matter is pressureless, w m = 0, the total 
energy density obeys the conservation equation 



dpt. 



+ -(1 + LU^Sl^ptot = . 



da a 

Therefore, assuming a constant w^, we obtain 

Ptot = Poa~ 3 [l - fi0,o(l - a 4 )]" 3 ^ , 
and the Friedman equation 



Hna 



5 [l-fVo(l-a«)]- 



(6) 



(7) 



(8) 



It is worth remarking that this expression for H 2 cannot 
be reproduced by a simple model of varying w, as e.g. 
w(z) = wq + w\z. This means that the effect of the 
coupling is intrinsically different from the effect of a time- 
dependent equation of state. 

The relation between the scaling behavior and the cou- 
pling S is easily derived by imposing the relation Q 



We find then 



t(-^) = °- 

dt p m ai 



(£ + 3w ) 



P,P + Pm 



(9) 



(10) 
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Substituting <@J we obtain the evolution of the cou- 
pling as a function of a 



So 



where 



(11) 



(12) 



Q<pfi + (1 - ^0,o)a 5 ' 

So = — f20,o(£ + 3w</>) . 
III. ADDING THE BARYONS 



So far we assumed a single matter component. How- 
ever, as already anticipated, if dark energy interacts with 
baryons then a new long-range force arises, on which 
there are strong experimental upper limits @,0j]|]. The 
baryonic component has therefore to be assumed ex- 
tremely weakly coupled; for simplicity, we assume here 
that the baryons are totally uncoupled. In principle, of 
course, we should allow the possibility that also part of 
the dark matter itself is uncoupled, as in Ref. but for 
simplicity we restrict ourselves to the basic case in which 
all dark matter couples with the same strength. Here 
we derive the corresponding formulae when uncoupled 
baryons are added to the cosmic fluid. As it will ap- 
pear clear, adding a small percentage of baryons at the 
present does not change qualitatively the fit to the su- 
pernovae; however, it changes dramatically the past and 
future asymptotic behavior of the cosmological model. 
In all the plots and numerical results we always assume 
fl b fl = 0.05. The parameter A now becomes 

A _ ^0,0 



1 — fi^o — " 



'0,0 



from which 



S (l - O bj0 ) 



^0,0 + (1 — ^6,o — ^0,o)a * ' 
so that the present coupling is 



(13) 
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6.0 



-(£ + 3u; ) 



(14) 



As it appears from Fig. 1, when £ > the coupling varies 
from zero in the past to a constant value S± = —(^ + 3w^,) 
in the future. As expected, the standard model of 
uncoupled perfect fluid dark energy is recovered when 
£ = — 3w0. For £ < the behavior is opposite. On the 
other hand, models with a constant S (investigated in 
[l4T p are obtained for £ = 0. In this case, ~ p m ; such 
behavior has been called "tracking" dark energy when 
the regime is transient an d "stationary dark energy" 
0,01 when the regime is the final attractor. To empha- 
size the connection between £ and the coupling, we will 
use in most cases the present coupling So = 6 (a = 1) as 
free parameter instead of £. The parameters that char- 
acterize the model are therefore (Q(j>,o, w^, So)- ACDM 



corresponds to w^, = —1 and So = 0. The trend of the 
coupling 5 can be expanded for low redshifts as 

X X { ~\ £.^m,0 ■. 



6,0 



The sign of £ is therefore also associated to the sign of 
the time derivative of the coupling. 
The Friedman equation is now 



H 2 = HqE(o) 



where 



E 1 



(l-fi M )a- d [l 



0,0 



(1-a 



£vr 3 - 



1 — £\o 

which for £ = becomes 

H 2 = H 2 {{1 - n b ,o)a- p + fl bfi a- 3 } , 



(15) 

(16) 
(17) 



where (3 = Iz^M+^Ag _ 

' 1 — "6,0 

It is interesting to derive the asymptotic limits of 
^6,m.0 (Figs. 2-4) Notice that, assuming < 0, the 
asymptotic future and past behaviour of the total den- 
sity and of the density parameters of all the components 
depends only on the parameters £ . For £ > we have 
the limits for a — > 
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(1-O 6 , )[l 



1-Qb 



-3uv/€ + ft 



6,0 



(18) 

(19) 
(20) 



while for £ < we have Cl b — ► 1 and Cl m , f2^ — > 0, so that 
the baryons dominated the past evolution even if their 
present density is very low. The future asymptotics for 
£ > is instead always dominated by dark energy, 



(21) 
(22) 
(23) 



n b - 


+ o, 
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1, 



For £ < we have 



n h 



^6.0 



fi?T 



(l-0 M )[l- T %£-]-3^ + fi M 

0, 



and finally for £ = 0, 

fib -> 0, 

^ 1 



0,0 



1 — ^6,0 ' 
fi<4,0 



i - n 



6.0 



> (24) 

(25) 
(26) 



(27) 
(28) 

(29) 
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It could be assumed that the negative values of £ are to 
be excluded, since they imply absence of dark matter in 
the past. However, we are here studying the model only 
in a finite range of z near the present epoch, and for this 
reason even negative £ cannot be excluded a priori. More- 
over, baryons dominate in the past for £ < and they 
could conceivably drive fluctuation growth. However, we 
will see that in fact values £ < are not favoured by SN 
data. 

There are other interesting features to remark in the 
general behavior of the density fractions. First, for So ^ 
the ratio f2(,/f2 m varies in time. If £ > 0, the ratio 
varies from a constant value in the past to infinity (if 
So < 0) or zero (if So > 0) in the future. If £ < the 
ratio decreases always (for w<p < 0): this variation in the 
baryon-to-matter ratio due to the dark matter coupling 
could provide additional testing ground for the coupling 
(e.g. 0). 

Second, for £ < the dark energy density vanishes 
both in the past and in the future, for all values of the 
equation of state: this means that the acceleration is only 
a temporary episode in the universe history, as shown in 
Fig. 4. The universe was dominated by the baryons in 
the past, by dark energy at the present and by a mix of 
baryons and dark matter in the future (for the param- 
eters employed in Fig. 4 the final value of fib is very 
low and cannot be distinguinguished from zero). For in- 
stance, assuming = — 1, fi^ = 0.7, £ = — 1 it turns out 
that the acceleration ends in the future at z = —0.85. 
This is a particularly striking example of how the cou- 
pling might completely modify the past and future be- 
havior of the cosmic evolution. 

Another example comes from the existence of the fu- 
ture singularity known as "big rip", i.e. an infinite growth 
of the total energy density in a finite time ^l- From the 
Friedmann equation ill fit we see that if £ > 0, for a large 
scale factor H 2 behaves as a^ 3 ^ 1 "^, as in the standard 
uncoupled model, so there is a big rip if w<j, < —1. If, 



instead, £ < 0, then H 



. This means that a nega- 



tive £ prevents the big rip for any value of w^. From 11411 
we see that, assuming f^o -C 1, the big rip is prevented 
when 




Figure 1: Behavior of S(z) for fi^ = 0.7 and w,p = —1. Full 
line: £ = —4; dotted line: £ = 4; dashed line: £ = 0. 
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Figure 2: Behavior of Q(a) for £ = 4 and w$ = —1, fixing 
flb,o = O.05,f2m,o = 0.25 and Q.4, = 0.7. Full line: dark energy; 
dotted line: baryons; dashed line: dark matter. 



So > 3|w0|fi0,o 

However, we'll find that this region of parameters space 
is rather disfavoured by SN data. 

We are interested also in evaluating the epoch of ac- 
celeration. The acceleration begins when a = i.e. for ( 



* = i^ { (Q M -i)[i- J^ ( i-^)]- 3 ^-n M 

a z 1 — s it, 



^ t — "6,0 



6,0 

(l-a«)]- 3 ^~ 1 } = 



(30) 



which will be solved numerically later on. For £ = the 
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Figure 3: Behavior of fi(a) for £ = and = — 1 (the other 
parameters are as in the previous figure). Full line: dark 
energy; dotted line: baryons; dashed line: dark matter. 
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Figure 4: Behavior of S7(a) for £ = —4 and w,p = — 1 (the 
other parameters are as in the previous figure). Full line: 
dark energy; dotted line: baryons; dashed line: dark matter. 
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Figure 5: Behavior of Sl(a) for the best-fit values of £ and 
(the other parameters are as in the previous figure). Full line: 
dark energy; dotted line: baryons; dashed line: dark matter. 



solution is 

1 — f4,o — 3^0,^0 

Zacc — [ 7^ J ' J- • 

"6,0 

Before starting with the data analysis, we need still 
another equation, the age of the universe. This is given 
by 



1 f 1 da 



to aE (a) 
where E(a) is given in Eq. I|16ll . 

IV. COMPARING WITH SNIA 

We used the recent compilation of SNe la data 0] 
(gold sample) to put constraints on the parameters 



^m,o, w <t>i &o entering the expression for H(a). In all sub- 
sequent plots we marginalize over a constant offset of the 
apparent magnitude, so that all results turn out to be in- 
dependent of the present value of the Hubble parameter. 
In all the calculations we fix f^o = 0.05. 

The overall best fit is fi^o = 0.62, w^, = —1.9, Jo = 
-1.5 (see Fig. 5) with a x 2 =173.7 (for 157 SN), to be 
compared to a \ 2 = 178 for the flat ACDM. In Fig. 
6 we present the main result of this paper: the con- 
tours at 68%, 95% and 99% of the likelihood function 
on the parameter space So,Wj>. The third parameter, 



0,0, 
<!>.() 



has been marginalized over with a Gaussian prior 
= 0.7 ± 0.1. The ACDM model lies near la from 
the best fit. The line £ = is tangent to the 3a con- 
tour: negative values of £ lie above this line and appear 
to be very unlikely with respect to positive values. It is 
impressive to observe how the vast majority of the like- 
lihood lies below the uncoupled line So, which is almost 
tangential to the la contour, and leftward of w^, = — 1. 
Taken at face value, this implies that the likelihood of 
a negative coupling (dark matter mass decreasing with 
time) contains 99% of the total likelihood, and that the 
likelihood of a phantom dark energy (w$ < —1) is 95% of 
the total likelihood. Our conclusion is therefore that the 
SNIa data are much better fitted by a negatively coupled 
phantom matter than by uncoupled, positively coupled 
or non-phantom stuff. It is interesting to note that it has 
been recently observed [23] that coupled phantom energy 
induces a repulsive interaction (regardless of the sign of 
the coupling). 

In Figs. 7 and 8 we plot the 1-dimensional likelihood 
functions for So, w^,, marginalizing over the other param- 
eters. The main conclusion is that current supernovae 
data exclude a coupling So > at 99% c.l. if no priors 
are imposed on w$. The lower bound is much weaker: a 
negative So down to So = —7 is allowed to to 68%. At 95% 
level, the formal lower limit to So is —10 but the conver- 
gence to zero is very slow and values very large and nega- 
tive of Sq cannot be excluded. In terms of the parameter 
£ we obtain £ > at 95% c.l. while essentially no upper 
limit can be put with the same confidence. In terms of 
m'/m, this implies that the dark matter mass can vary 
by a fraction <5o in a Hubble time. Looking at Fig. 6 one 
sees that a zero or positive coupling is instead preferred 
if wo > — 1- Applying the prior w<f, > — 1 the limits on 
So narrow and move to higher values: —2.5 < So < 1.2 
(95%). 

From Fig. 8 we see that — 3 < wj, < — 1.2 at 68% and 
—4.2 < Wcf, < — 1 at 95% c.l.. In the same Fig. 8 we plot 
the likelihood for w^, in the uncoupled case. As it can 
be seen, the coupled case extends considerably the al- 
lowed region of w$ , especially towards large and negative 
values: a negative coupling favors phantom equation of 
states. This complements the results of ref. [lj| m which 
it was found that a positive coupling (corresponding to a 
stationary behavior with £, — and < 0) favours non- 
phantom matter. This is a clear-cut example of how the 
conclusions regarding the nature of dark energy depend 




Figure 6: Likelihood contours at 68%, 95% and 99% c.l., in- 
side to outside, marginalizing over fl^ = 0.7 ±0.05 (Gaussian 
prior). In this and the following two-dimensional plots the 
star marks the best fit; the horizontal line indicates the de- 
coupled models, the vertical line separates the phantom (on 
the left) from the non-phantom models (on the right); the 
ACDM model is at the crossing point of the two lines. Below 
the dashed line £ > (fixing f2</, = 0.7). On the axes, the grid 
we used for the computation. 



Figure 7: Full line: likelihood for So (marginalized over 
and w$). The horizontal dashed lines give the 68%, 95% and 
99% c.l., top to bottom. Dot-dashed line: likelihood for So 
with the prior > — 1. 



crucially on its coupling to the rest of the world. It is in- 
triguing that the correlation between So and w^, crosses 
approximatively the ACDM case So = 7 w<f, = —1: al- 
though not particularly favoured, ACDM remains a per- 
fectly acceptable model also with respect to coupled dark 
energy. 

We now impose on the 2D likelihood the contour lev- 
els of the acceleration ijiffljl . Since a negative So means 
a more recent surge of the dark energy it is to be ex- 
pected that the likelihood favours a recent acceleration. 
It turns out that indeed the best fit corresponds to an ac- 
celeration epoch z acc ~ 0.3; however, at the same time, 
an acceleration z acc > 1 lies near 2a and cannot be ex- 
cluded with large confidence (Fig. 9), especially if one 
excludes phantom states of matter. 

The age contours are compared to the likelihood in 
Fig. 10. Most of the likelihood lies within the acceptable 
range 11 and 14 Gyr (for h = 0.7 and Q^q = 0.7). The 
age constraints are therefore rather weak. 

We can also marginalize on So and plot the likelihood 
for ttmOjWcf, (see Fig. 11). In Fig. 12 we compare the 
likelihood for Q m0 for the uncoupled case, So = and 
for the general case. Again the result is that the likeli- 
hood for ttmo widens considerably. Now practically any 
value from fl m o = to f2 m o = 0.7 is acceptable, with 
a broad peak around 0.25, while in the uncoupled case 
f2 TO o peaks rather tightly around 0.4 (remember that we 
are marginalizing over all values of w^); if we restrict to 
Wcf, > —1 then we have roughly 0.1 < Q m o < 0.3 at 95%. 
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Figure 8: Full line: likelihood for w$ (marginalized over fi^, 
and So). The horizontal dashed lines give the 68%, 95% and 
99% c.l., top to bottom. Dotted line: fixing do = 0. 



V. FORECASTS FOR FUTURE ESPERIMENTS 
A. SNAP 

There are several project to extend the SNIa dataset 
both in size and in depth. No doubt they will increase 
our understanding of the dark energy problem. Here we 
investigate the potential to constrain the parameters of 
the interacting model with a SNIa dataset that matches 
the expectation from the satellite project SNAP. 

We generate a random catalog of redshifts and magni- 
tudes of 2000 supernovae from z = to z = 1.7, with a 
r.m.s. magnitude error Am = 0.25, distributed around a 
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Figure 10: As in Fig. 6 with overimposed the lines of equal 
age T = 11, 12, 14, WGyr (assuming VL^ = 0.7, h = 0.7) , 
bottom to top. 



Figure 12: Likelihood for Q, m marginalized over 5o,w$. The 
dotted curve is for uncoupled dark energy (So = 0). 



ACDM cosmology with Slmo = 0.3, So — 0, wq — — 1. The 
redshift have been assumed uniform in the z-range; very 
likely this is not a good approximation to what SNAP will 
produce but it is rather difficult at the present stage to 
predict what the final z-distribution will be. It is likely in 
fact that the statistical significance of the larger number 
of expected sources at high z will be at least partially 
reduced by additional uncertainties like lensing effects, 
redshift errors, spread in the calibration curve etc.. So 
we preferred to keep the forecast as simple as possible 
in order not to introduce additional, and not well moti- 
vated, parameters. 

In Fig. 13 we show the likelihood marginalized over 
f2 m o- The 68% errors are of order 0.2 for w and 0.5 for 
6q . Such an experiment will therefore be able to increase 
the precision in wq in So by a factor of five, roughly. 



B. Baryon oscillations 

The main limit of the SNIa method is that even satel- 
lite experiments like SNAP do not expect to detect a sig- 
nificant number of sources beyond z = 1.7. It has been 
suggested that an interesting possibility to probe deeper 
the cosmic history is the reconstruction of the baryon os- 
cillations in the galaxy power spectrum |2lJ ■ The method 
exploits the wiggles in the power spectrum induced by the 
acoustic oscillations in the baryon-photon plasma before 
decoupling as a standard ruler. When viewed at different 
z's, the size of the oscillations map into angular-diameter 
distances that probe the cosmic geometry just as the lu- 
minosity distances of the SNIa. 

Here we forecast the constraints on , S by baryon 
oscillations assuming exactly the same experimental 
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ble I we list all present and future contraints (at 68%) we 
derived in this paper, bearing in mind that the future ex- 
periments have been tested against a ACDM target only: 
a different target cosmology implies in general different 
errors. 
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AS 


Riess et al. (2004) SNIa 


ps 1 


ps 2.5 


0scillations:2 = 1+2 = 3, spectr. 


0.22 


0.09 


Oscillations: 2 = 1+2 = 3, phot. 


0.57 


0.54 


Oscillations:2 = 1, spectr. 


0.51 


0.81 


0scillations:2 = 1, phot 


0.78 


2.52 


Oscillations:2 = 3, spectr. 


0.36 


0.10 


Oscillations: 2 = 3, phot. 


0.64 


0.66 


SNAP 


ps 0.2 


ps 0.5 



Table I 



Figure 13: Likelihood contour for a SNAP-like experiment, 
assuming ACDM (fl m — 0.3, So = 0, = —1) as target 
cosmology. 

specification of ref. j^l, to which we refer for all the 
details. To summarize, we evaluated the Fisher ma- 
trix of several combined datasets: five deep surveys of 
200 deg 2 each binned in redshift and centered at z = 
0.7, 0.9, 1.1, 1.3, 3, plus a survey similar to the Sloan Dig- 
ital Sky Survey (SDSS) in the range 2 = — 0.3, plus a 
cosmic microwave (CMB) experiment similar to the ex- 
pected performance of the Planck satellite. The five deep 
surveys are assumed either spectroscopic (negligible er- 
ror on z) or photometric (absolute error Az = 0.04). 
The reference cosmology is again ACDM as above. With 
respect to the method of [23 we marginalize over the 
growth function. 

The results are shown in Fig. 14 for various combina- 
tions. The label 2 = 1 denotes the combination of the 
four surveys at z = 0.7 — 1.3; 2 = 3 denotes the furthest 
survey. In all cases we include the Fisher matrices of 
SDSS and CMB. As expected, the deepest data are the 
most powerful: the dipendence of the angular diameter 
distance D on So (i.e. <91ogZ?/<91og(5o) increases by more 
than a factor of 2 from z = 1 to z = 3. Consider first 
the spectroscopic case (continuous curves): including the 
survey at z = 3 pushes down the 68% c.l. errors on 
to roughly 0.36 (not far from SNAP forecasts) and to 0.1 
the error on So, quite better than the SNAP forecasts. 
Adding the surveys at z ps 1 almost halves the error on 
w<f, (Awq ps 0.22) while is not particular effective versus 
the error on Sq. On the other hand, the photometric sur- 
veys give constraints Aw,/, ps 0.6, ASo ~ 0.5. Overall, 
we conclude that the baryon oscillation method at z ps 3 
improves upon the SNAP experiment for as concerns So, 
while has a similar efficiency in constraining w<j,. In Ta- 



VI. CONCLUSIONS 

This paper studied the behavior of a model with three 
components, baryons, dark matter and dark energy. Con- 
trary to most similar analyses available in literature we 
included a general phenomenological coupling between 
dark matter and dark energy specified by a simple but 
rather general scaling relation. By this extension of 
dark energy models we have been able to put bounds 
on the present interaction between the two components 
or, equivalently, to the rate of change of the dark matter 
particle mass. Moreover, we can appreciate the extent to 
which the constraints on such fundamental quantitites as 
fl m and Wj, depend on the assumptions concerning the 
dark energy interactions. 

We obtained several results that we summarize here. 

1. We find the constraints on the coupling constant 
So, which can be interpreted as the rate of change 
of dark matter mass per Hubble time. We obtain 
-10 < S < -1 at 95% c.l. and S < at 99%. 

2. We find -4.2 < w < -1 (95% c.l.). Together with 
the bounds on So we conclude that the current SN 
data set favours negatively coupled phantom dark 
energy. In general, we find that the equation of 
state correlates with the coupling: positive cou- 
pling implies w^, < —1, negative coupling implies 
phantom energy. Quite remarkably, imposing zero 
coupling peaks the cosmological constant as pre- 
ferred value of u>0 . 

3. Marginalizing over the coupling we derive new 
bounds on £l m o- In particular, we find —4 < w^, < 
-1 and 0.05 < f2 m0 < 0.65 (95%) . The dark 
matter density allowed region widens considerably 
when a nonzero coupling is introduced in the model. 
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4. We find as best fit a value for the beginning of ac- 
celeration z acc = 0.3 but we cannot exclude earlier 
acceleration (z > 1) at more than 95% c.L. 

5. We find that the addition of even a small fraction 
of uncoupled matter, e.g. baryons, modifies pro- 
foundly the asymptotic behavior of the model. In 
particolar, for £ < we find models in which the 
epoch of dark energy dominance (and acceleration) 
is only temporary. Notice that this occurs for w^, 
negative and constant. 

6. The ratio of baryons-to-dark matter varies in our 
model. Depending on £ it may increse or decrease 
with time. This may offer new methods of con- 
straining a preferential coupling to dark matter. 

7. The big rip that in uncoupled models occurs for 
Wtj, < — 1 can be prevented if So > 3\wtj,\fl^fll how- 
ever, these values are unlikely at more than 99% 
c.L 

Future experiments will constrain the equation of state 



and the coupling to a much better precision. We find 
that an experiment with the specification of SNAP might 
reduce the errors on So,w<p by a factor of five roughly. 
The method of the baryon oscillations could reduce the 
error on So by a factor of 25 roughly. Needless to say, 
these forecasts depend on the exact experimental setting; 
however, they give a feeling of the expected precision on 
the dark energy parameters that can be reached a few 
years from now. Whether the final outcome will show 
any trace of coupling (or, for that matter, any trace of 
deviation from a pure cosmological constant) is one of 
the most exciting question to ask. 
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Figure 14: Likelihood contours at 68% estimated from a 
Fisher matrix analysis of baryon oscillations of several com- 
bined surveys. In all plots the inner continuous curve is for 
spectroscopic surveys, the outer dashed curve for photometric 
surveys with redshift error 0.04. Top: four survey at z « 1; 
middle: survey at z ~ 3; bottom: combined surveys at z w 1 
and x^3. In all cases, the SDSS survey and the Fisher matrix 
for a CMB experiment similar to Planck has been included. 



